Controlling the dynamics of an open many-body quantum system with localized dissipation 
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We experimentally investigate the action of a localized dissipative potential on a macroscopic matter wave, 
which we implement by shining an electron beam on an atomic Bose-Einstein condensate (BEC). We measure 
the losses induced by the dissipative potential as a function of the dissipation strength observing a paradoxical 
behavior when the strength of the dissipation exceeds a critical limit: for an increase of the dissipation rate the 
number of atoms lost from the BEC becomes lower. We repeat the experiment for different parameters of the 
electron beam and we compare our results with a simple theoretical model, finding excellent agreement. By 
monitoring the dynamics induced by the dissipative defect we identify the mechanisms which are responsible 
for the observed paradoxical behavior. We finally demonstrate the link between our dissipative dynamics and 
the measurement of the density distribution of the BEC allowing for a generalized definition of the Zeno effect. 
Due to the high degree of control on every parameter, our system is a promising candidate for the engineering 
of fully governable open quantum systems. 

PACS numbers: 67.85.Hj, 03.75.Gg, 34.80.Dp 



Gathering information from a quantum system is never free 
of cost. Every measurement process provides a coupling be- 
tween the quantum system and the (classical) environment, 
which leads to non-unitary dynamics, and in some cases to the 
destruction of essentially quantum effects. The elusive transi- 
tion from the quantum to the classical realm must therefore be 
inherent in the processes that the environment induces on the 
system. In recent decades several advances have been made 
in the study of environmentally induced phenomena like deco- 
herence and decoherence-induced selection of preferred states 
(einselection) lUrS . More recently environmental action has 
been used to manipulate qubits in a system of trapped ions fl. 
The knowledge and the mastering of the action of the envi- 
ronment are essential for taming errors in quantum computa- 
tion schemes or to engineer decoherence-free subspaces 
for qubits [7-9], and are also key to understanding the emer- 
gence of classicality from the quantum |2l|3|]. In the context 
of the theory of open quantum systems, environmental action 
gives rise to effective Hamiltonians which can contain imag- 
inary terms Since these terms actually arise from 
a collection of an enormous number of degrees of freedom 
|3l], however, a complete experimental control over them ap- 
pears overly challenging. Here, we report the engineering of 
a fully controllable, environmentally induced imaginary po- 
tential acting on a quantum system, and present observations 
of the subsequent induced dynamics. The localized imagi- 
nary potential is realized by the almost pure dissipative action 
of an electron beam (EB) on an atomic Bose-Einstein con- 
densate (BEC). We show that such a potential can be used to 
describe a continuous measurement process that can exhibit 
a generalized version of the so-called Zeno effect. The com- 
bination of the robust and macroscopic many-body quantum 
behavior of a BEC and the high tunability and precision of 
the EB promotes such a system as a paradigm for governable 
open quantum systems. 

One of the most striking properties of BECs is that, de- 
spite their many-body nature, they can be described to a good 



approximation by a mean-field wavefunction obeying the so- 
called Gross-Pitaevskii equation (GPE). This remains valid 
also when the BEC is coupled with the environment. Starting 
from the Lindblad master equation ihd,p — [H,p] + ihJ£p, 
where p is the density operator of the many-body system, H 
is the Hamiltonian operator, and Jz? is the dissipation operator 
such that Jz^p = - j't/xy(x)/2[ v P+ l Pp + p*P +l P - 2 l Pp v P + ], 
with y(x) the local dissipation rate, we can write the equation 
of motion for the expectation value of the bosonic field opera- 
tor 'P as dr^P) = TrQVdtp), which leads to a time-dependent 
GPE with an additional imaginary term ifllUl^l : 



dy(x,f) 
dt 



h 2 v 2 

2m 



-V ext +g\y{x,t)\ 2 -ih 



\j/(x,t). 

(i) 

Here y/ is the BEC wavefunction, obeying the constraint 
/ 1 y/(x,?)| 2 c/x = N(t), V ex t the trapping potential and g = 
Anh 2 a/m, a being the s-wave scattering length. Notably our 
technique allows independent control of the Hamiltonian and 
dissipative terms of this equation. The ability to describe our 
open quantum many-body system with such a simple expres- 
sion is a key asset for understanding and mastering its dynam- 
ics. 

Experimental implementation In our experiment we pre- 
pare a pure BEC of 75x1 3 atoms in a single-beam optical 
trap by means of forced evaporation. Once the evaporation 
is over, we shine a focussed EB right at the center of the 
BEC. The EB is produced by a commercial electron micro- 
scope mounted inside the vacuum chamber 01411 . The electron 
microscope is able to generate a beam of 6 keV electrons with 
variable beam extensions and currents. When the electrons 
impact on the BEC they collide locally with the atoms, ion- 
izing or exciting them. Those atoms which have undergone 
an electron collision escape from the trapping potential (Fig. 
la). The EB thus locally dissipates the BEC. The ionized 
atoms, roughly 40% of all those scattered, are then directed 
to an ion detector, where their arrival times are registered. 
While escaping from the trapping region the ions can col- 
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FIG. 1 : (Color online) a) The electrons locally collide with the atoms 
constantly dissipating the BEC. b) Temporal resolved signal from 
the ion detector. The bin size is 1 ,us. Points are experimental data 
averaged over 1800 experimental repetitions, while the solid curve is 
the numerical simulation (see text). After 5 ms we typically collect 
~ 450 ions. 

lide with the trapped atoms producing additional losses lUUl . 
The total detection efficiency 77 is the product of the branch- 
ing ratio (40%) and of the combined ion optics and detector 
efficiency (75%). Details of the experimental apparatus can 
be found in 11141 1 1 511 . If the EB is rapidly moved in a con- 
trolled pattern, the whole column density profile of the BEC 
can be reconstructed lfl4lll6ll . Here we keep the EB fixed in 
the center of the BEC, and monitor the subsequent induced 
dynamics by looking at the temporal signal from the ion de- 
tector (Fig. lb). By controlling the beam parameters, we can 
engineer the dissipative term in equation (HJ: we write it as 
y(x) (2new 2 ) exp(- (x 2 +y 2 ) /2w 2 ), I being the EB cur- 

rent, <7 the electron-BEC scattering cross-section [11], e the 
elementary charge and w the standard deviation of the spatial 
electron distribution, assumed to be gaussian 

Comparison between experimental results and theoretical 
expectations In Fig. 2 we report the number of ions col- 
lected in the first 5 ms of continuous dissipation as a function 
of the EB current for three different values of w. Notably we 
observe that the number of ions produced, as a function of 
the EB current (i.e., of the number of electrons sent on the 
atoms), shows a non-monotonic dependence. In other words, 
starting from a critical value of the EB current, the harder we 
try to dissipate, the less we manage to do it. This paradoxi- 
cal behavior is more marked for smaller values of w. In the 
same figure the data are compared with the results obtained 
by numerically solving equation (Q]i, additionally taking into 
account secondary effects like ion-atom collisions 111 111 . The 
agreement is very good (and the same agreement is visible in 
Fig. lb), demonstrating that the description of the EB as a 
pure dissipative potential is sufficient to capture the observed 
main features. A detailed description of the dynamics that 
leads to the curves reported in Fig. 2 will be given in the fol- 
lowing. From simple textbook calculations, or from more for- 
mal analysis like the one made in 11711 . it is easy to verify that a 
localized imaginary potential U induces total reflection as the 
strength of the potential goes to infinity. Hence the effective 
quantum dissipation vanishes when the localized imaginary 
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FIG. 2: (Color online) Number of ions collected within the first 5 ms 
of continuous dissipation on a BEC as a function of the EB current 
/. The three panels report the data obtained with w = 127(5), 170(7) 
and 212 (8) nm, from left to right. Each data point is the average over 
75 experimental repetitions. The error bars are mainly due to shot- 
to-shot fluctuations in the overall ion detection efficiency. The solid 
lines are the number of dissipated atoms resulting from the numerical 
simulations (see text). Please note that in the sum also the initial 
decay visible in Fig. lb is included. The scale on the top reports the 
strength of the imaginary potential U /h = y(0)/2 corresponding to 
the measured current. 



potential is either zero or infinity, implying the existence of a 
maximum of dissipation for some finite value of U. This ex- 
plains on a qualitative basis the observed non-monotonicity. 
The position of the maximum is of special importance, since 
it sets the parameters which allow one to engineer the most 
efficient possible absorbing potential. As an example, in lfl7ll . 
where the time of arrival of a Id wave packet is measured by 
a step-like potential, the maximum dissipation is analytically 
calculated to be Um — 10. 6E, where E is the energy of the 
wave packet. In our case the presence of the non-linearity 
and the less idealized conditions do not allow for an ana- 
lytic solution, but from the numerical results plotted in Fig. 
3a we have found Um — 8/xexp(— w/d), where and d are 
respectively the chemical potential and the healing length of 
the unperturbed BEC. From Fig. 3a it also appears that in- 
creasing the size of the EB not only moves the position of 
the maximum dissipation to higher values of 7, but also in- 
creases the number of produced ions for a given current, and 
reduces or washes out the effect of the reflection. Clearly, 
when U = ny(0)/2 > Um, a decrease of the probe size (an in- 
crease of the resolution) leads to a lower production of ions, 
making the system more resistant to the environmental action. 

Comparison to classical systems To ascertain to which 
extent our observations are peculiar to the wave nature of the 
BEC, we have repeated the experiment on a thermal gas of 
4x 10 s atoms at 1/iK with a beam of w = 170 nm. The results 
are reported in Fig. 3b, where a simple monotonic behav- 
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FIG. 3: (Color online) a) Theoretical curves of the number of ions 
produced in 5 ms as a function of U/h solving equation l[T) for dif- 
ferent values of w. The values of Um obtained using the approximate 
expression given in the text are shown as open diamonds over the 
corresponding curves, b) Number of ions measured after 5 ms of 
dissipation for a thermal cloud as a function of the EB current, with 
w — 170(7) nm. The solid line is the result of the corresponding nu- 
merical simulation using the molecular dynamics method, c) Com- 
parison between the theoretical curves of the number of produced 
ions as a function of / for the BEC and the corresponding (see text) 
classical analogue (w = 170(7) nm). 



ior is observed. Such data are well reproduced by a classical 
molecular dynamics simulation, which includes the dissipa- 
tion induced by the EB [11]. We then extend the classical 
simulation to an atomic cloud which has the same density, 
number of atoms and trapping frequency as our BEC. Even 
though this does not represent any real physical system, it 
is instructive to compare the behavior of a quantum system 
with its hypothetical classical analogue. This comparison is 
made in Fig. 3c, where the monotonicity of the classical case 
is confirmed. From this we can conclude that the observed 
non-monotonicity is a purely quantum effect stemming from 
the macroscopic wave nature of the BEC II32I1 . Moreover it 
is evident that the effect of the quantum reflection from the 
imaginary potential leads to a suppression of dissipation in the 
quantum case, which is already notable for very low currents. 

Dissipative dynamics In order to gain a deeper insight 
into the dissipation-induced dynamics, we now look in de- 
tail at the time-resolved signals coming from the ion detector, 
reported in Fig. 4a. Initially the number of ions produced is 
well described by the exponential decay exp(— t y), /being the 
effective dissipation rate 111 ill . In Fig. 4b we show the integral 
of the signals in the first 5 fls, together with the simulated val- 
ues as a function of the EB current. In this phase, where no 
paradoxical behavior is either observed or expected, the EB 
bums a hole in the BEC wavefunction (see Fig. 4d), defining 
a clear border between the space "inside" the hole and that 
"outside". Thereafter, the number of ions produced becomes 
almost constant, signaling the onset of a quasi-stationary dy- 
namics 111 8fl . In this second phase, the reaction of the quan- 



tum system to the external perturbation takes place. When the 
strength of the dissipation is increased, the "outside" wave- 
function passes from a situation of almost total transmission 
to a situation where reflection takes the leading role. The 
non-monotonic dependence on the dissipation strength then 
becomes apparent (Fig. 4c). This represents the first exper- 
imental observation of the so-called back-flow paradox II 1911 . 
i.e., of the onset of a temporary reflection from a localized 
perfect absorber. The curves plotted in Fig. 2 and 3 are then 
the sum of different contributions like those in Fig. 4b and 4c. 

Dissipation as continuous measurement Finally we 
demonstrate that the controlled dissipation is equivalent to a 
local measurement of the BEC density, i.e., of the squared 
modulus of its wavefunction. Starting from equation ([TJ and 
defining (j> = yz/VN, where N is the number of atoms, after 
some algebra we obtain the equation 



dN{t) 
dt 



= -N(t) / y(x)\<i>(x,t)\ l dx. 



(2) 



The number of ions produced in a time interval At around a 
certain time t is ANj(t) = TJ J^^j^ \dN(t)/dt\dt. Hence we 
can conclude that what we perform is a direct measurement 
of the BEC density 1^(0 1 2 m me region illuminated by the 
EB, as in ffl [lrjll . Since the seminal trilogy on the time of 
arrival in quantum mechanics [17], imaginary potentials have 
been linked to the action of a measurement apparatus while 
later refinements [20] formally demonstrated the equivalence 
between a pulsed measurement with period 8t and a contin- 
uous dissipative potential U, provided that 8t ~ h/U. Our 
findings represent the experimental verification of the equiv- 
alence between the action of an imaginary potential and the 
one of a measurement apparatus on a quantum system. In- 
deed we show that a continuous measurement of the BEC 
density is nicely reproduced by introducing an imaginary po- 
tential in the corresponding Schrodinger equation. Further- 
more it is well understood theoretically [|21 22], and verified 
experimentally 11231 - 12511 . that performing continuous measure- 
ments strongly modifies the pre-existing dynamics of a quan- 
tum system. This is known as the Zeno (or inverse-Zeno) ef- 
fect. In general, measurements are mainly performed on non- 
decaying systems and an extension of the standard definition 
to such systems is needed. In our case no pre-existing dynam- 
ics is present, since in the absence of the EB, the BEC is at 
rest. We have shown that the action of the continuous dissi- 
pative potential, or of the continuous measurement, strongly 
modifies the output of the measurement itself. In analogy with 
the standard definition we define dissipation induced Zeno dy- 
namics (DZD) when dNi(t)/dU < 0, N{(t) being the number 
of ions produced in the time t , and dissipation induced simple 
dynamics (DSD) when dNi(t)/dU > 0. These definitions ap- 
pear to be the natural extension of the standard ones, since the 
onset of the DZD requires large values of U = hy(0) /2, which 
corresponds to pulsed measurements with small 8t [20]. From 
the definitions it follows that the DSD is observed where the 
dynamics is dominated by the Hamiltonian term of the Lind- 
blad master equation, while the onset of the DZD corresponds 
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FIG. 4: (Color online) a) Temporal resolved signal of the arrival time of the ions on the ion detector for different values of the EB current / for 
w = 170(7) nm. In order to enhance readability the data have been plotted with a binning of 50 /is. The insets b) and c) show the integrals of 
the signal in the shaded areas (the first 5 /is and from 1 to 5 ms respectively) for different values of / together with the theoretical calculations 
obtained solving equation (0. d) Points: scanning electron microscopy image of the BEC profile along the weak confining axis of the optical 
trap. The profile is the integrated column density along the direction of propagation of the EB. The scan is made after 1 ms of dissipation. The 
EB parameters are / = 150 nA and w = 106(5) nm. The depletion of the density is visible in the origin, i.e., in the center of the BEC. The solid 
line is the profile obtained numerically solving equation (0. 



to a dynamics governed by the dissipative term. We note that 
an effect resembling the DZD has been observed also in a 
system of decaying molecules in Id l26ll and in an attractive 
Mott-insulator state ll27ll . 

Conclusion and Outlook We have experimentally demon- 
strated the implementation of an open many-body quantum 
system whose Hamiltonian and dissipative dynamics can be 
independently and accurately controlled. In the case of ex- 
tremely strong and localized dissipation this can lead to the 
creation of dissipation-resistant states. The possibility to cre- 
ate such states in a controlled fashion can give new insights 
for engineering generalized environmental dark states. These 
kind of states are of fundamental interest and can possibly 
have practical applications in quantum computation schemes 
112811 . And inasmuch as our technique exploits the demon- 
strated link between dissipation and measurement, it can be 
used to address fundamental issues in quantum mechanics, 
like the definition of the time of arrival [17]. The dissipation 
mechanism studied in the present paper is also particularly 
suited for lattice systems [29], thanks to its localized charac- 
ter and hence to the ability to selectively control the dissipa- 
tion in a single lattice site. Indeed the use of the EB offers the 
unique possibility to create and study long-living exotic states 
in optical lattices II 1 3fl and to characterize the interplay be- 
tween dissipation and interactions [ 11 , 3(J, and so would give 
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SUPPLEMENTARY MATERIALS 

From the Lindblad Master Equation to the Gross-Pitaevskii 
equation We begin with a master equation of the Lindblad 
formJl 



p l P t (x)4 , (x) - 24 , (x)p4 ,+ (x)) . 



(5) 



iti'p = [H,p] +ihJfp, 



(3) 



which describes the time evolution of the density operator of 
the system under the combined action of the Hamiltonian 

H = Jd\v\x) ^|-V 2 + V ext (x)^ + (4) 

|y c/ 3 jc4 ,t (x)4 ,t (x)4 , (x) , l'(x) 
and a local loss term, described by the dissipation operator 
i?p = - J rf 3 Jt^p ( l P t (x) l P(x)p + 



Here, ^(x) and ^(x) are the bosonic field operators, obeying 
the bosonic commutation relation 



[«P(x),«pt(x')] =5 3 (x-x'), 



(6) 



the interaction strength is denoted with g and V e xt(x) is an 
external potential. The loss term y(x) describes the rate of 
particle losses at the position x. Similar master equations for 
optical lattices with leaky sites have been studied in 11131, 
The time evolution of an operator is then given by the differ- 
ential equation 

0?{x))=Tr (7) 
Plugging equation (3) into (7) yields after some algebra 

h 1 



(8) 



s4rt(*)*(x)]*(x)>-^ <❖(*)>. 
Performing the usual mean-field approximation 



(9) 



we arrive at equation (1) of the main text, i.e., to the time- 
dependent Gross-Pitaevski equation with imaginary potential 
(IGPE). The comparison between the experiment and the nu- 
merical solution of the IGPE suggests that the mean-field ap- 
proximation is well justified for our parameters. 

Beam characterization We experimentally determine the 
value of (7 by monitoring the short-time decay of the ion sig- 
nal for a fixed current and different electron beam sizes. In- 
deed analysing the solutions of the IGPE we have found that 
such a decay is approximately exponential with a time con- 
stant T = Iff — l6ew 2 /Io. The EB size w is derived in two 
different ways. We first measure w by fitting with an inte- 
grated gaussian function the electron signal originating from 
the scan of the beam through the edge of an auxiliary metallic 
target. The position of the beam focus and hence the value of 
w at a fixed point can be easily varied by changing the cur- 
rent in its final focussing stage. A calibration of the beam size 
as a function of this current is also used to recover the mea- 
sured w when the working position of the microscope is finely 
adjusted at the actual location of the atoms. For comparison 
we also extract w directly using the atomic sample as a target, 
by measuring the visibility of a two-dimensional deep optical 
lattice imaged by the EB itself. The lattice that we employed 
has a spacing of 547 nm and a depth of 21 E r . The beam size 
w is obtained by deconvolving the density distribution of the 
atoms in the lattice by their expected Wannier functions distri- 
bution. Note that the EB size does not vary appreciably along 
the radial extension of the cloud (~ 7 fxra), as its Rayleigh 
length is always conveniently larger (> 30 fim). The two mea- 
surements yield compatible results and are also in accordance 
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can find that the EB generates a magnetic field which has the 
following expression: 



B e (R) 



HoIR -£ 

- - e zw 2 

Anw 1 



(11) 



The electron microscope machine produces itself a back- 
ground magnetic field Bf, g of ~ 1 G, directed along the di- 
rection of propagation of the EB and is always present at the 
atoms position. Due to this, the total magnetic field acting on 

the atoms is B(x) = ^JB 2 +Bj }g . Since B e is on the order of 

just 1-2 mG, it results that B ~ Bb g , making any effect due to 
the magnetic field produced by the EB B e negligible. 

We now consider the electric field induced on the atoms by 
the flux of electrons, which can be approximated by: 



FIG. SI: (Color online) Inverse decay time 1/t as a function of the 
EB current / for two different values of the beam size w = 170(7) 
nm (gray triangles) and w = 255(10) nm (red circles). The solid 
lines are linear fits to the data. The current / is measured by means 
of a Faraday cup, which is integrated in the experimental apparatus. 



with the measurement of 1/t as a function of the current / 
performed for two different beam sizes: w — 170(7) nm and 
w = 255(10) nm, see Fig. SI. From this analysis we obtain 
the value of the electron-BEC cross-section used in this work: 
a = 3.6(0.1) x 10~ 21 m 2 . 

Secondary Effects The main spurious effect is caused by 
the losses due to the secondary collisions of the produced ions 
with the trapped atoms. Indeed such ions, produced in the 
center of the BEC, travelling through the atomic cloud can 
collide with other atoms expelling them from the trap. This 
effect can in principle be reduced applying a strong electric 
field which extracts the ions faster; but it cannot be completely 
avoided. In order to take into account this loss channel the 
IGPE must be modified adding an ion-atom collision term in 
the imaginary potential such that: 



E e (R) 



2keqR\\ 



(1 



(12) 



hy(x) hy(x) ha(x) dNi(t) 
~~2 > 2 +_ 2 d~T 



(10) 



where dNi(t)/dt is the ions production rate and oc(x) = 
Aexp{-(x 2 +y 2 )/2W 2 ) is a gaussian function whose ampli- 
tude A and extension W depend on the extraction geometry 
of the ions. We have determined these two parameters fit- 
ting the time resolved curves like the one in Fig. lb, finding 
A = 1 x 10~ 3 and W = 20w. Taking into account this effect 
is important in our case mainly to estimate our ion detection 
efficiency, that is ~ 75%, but it has no significant impact on 
the physics presented. In the main text when we compare the 
theory with data we have included this effect while pure theo- 
retical graphs report the ideal case. 

Additional spurious effects that, in principle, must be taken 
into account are related to the real potentials stemming from 
the EB. Indicating with R the radial distance from the center 
of the EB in the direction perpendicular to its propagation, one 



where v e is the mean velocity of the electrons. This field leads 
to an effective potential orientated along the direction per- 
pendicular to the beam propagation Ue(R) — — l/2(XoE 2 (R) 
where ao is the ground state polarizability of the 87 Rb atoms. 
For the beam maximum current we employed in the present 
experiments, / = 500 nA, the ratio between the maximum 
modulus of the dissipative and electric potentials is ~ 580 and 
this latter contribution is again negligible. Finally, we have 
also taken into account for the effects of the real potentials of 
the EB in the numerical simulations finding, as expected, no 
appreciable correction to the calculations performed without. 

Molecular Dynamics With Dissipation The numerical 
simulations made to reproduce the data for the thermal cloud 
were performed using a molecular dynamics method. We have 
used a Runge-Kutta algorithm to simulate the dynamics of 
4x10 s atoms at the temperature of 1 fiK in a harmonic po- 
tential. For this simulation the atoms were supposed to be 
point-like, neglecting any kind of quantum effect like the fi- 
nite extension due to the deBroglie wavelength. Moreover we 
did not take into account the interactions between the atoms. 
The action of the electron beam was simulated introducing 
I/eAt electrons per m 2 , gaussian distributed in space along 
x and y, and uniformly distributed along z, the direction of 
propagation of the EB. At is the time step of the simulation. 
Every time the distance r between one atom and one electron 
is smaller than yj a j% the atom is supposed to collide with the 
electron leading to the loss of such an atom from the trap. We 
have found that this method, despite its simplicity, reproduces 
quite well our data except for a multiplicative factor, which 
takes into account the efficiency of the ion detector and the 
branching-ratio of the collision processes. 

Effect of the interactions We point out that a big advan- 
tage originating from the many-body nature of the BEC is the 
possibility to perform continuous measurements on a single 
wavefunction, giving direct access to the dynamics without 
the need of ensemble averages. At the same time the inter- 
actions between the atoms play a crucial role in determining 
the properties of such a wavefunction, through the non-linear 
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FIG. S2: (Color online) Calculated total number of ions produced in 
the first 5 ms as a function of the EB current for different values of 
the scattering length a. The position of the maximum dissipation is 
shifted to higher values of / as the interactions between the atoms are 
increased. The normalization is made to compensate the effect of the 
changing density (oc a~ 3 / 5 ). Clearly increasing interactions hinder 
the onset and the effects of the DZD. 



term in the IGPE, and have a direct impact on the dissipative 
dynamics. Even if experimentally we are not able to change 
the interactions we can do it when solving the IGPE changing 
the value of a. The results are reported in Fig. S2. 



As expecting a change in the interaction strength changes 
the position of the maximum of dissipation, since it depends 
on fi: the stronger is the repulsion between the atoms the more 
the onset of the DZD is shifted to larger values of / and the 
more its effects are decreased. Moreover an increasing in a 
implies a reduced production of ions, for a given y, due to the 
subsequent change in the density. As shown in Fig. S2, also 
a non-interacting BEC displays the non monotonic behaviour 
observed in the experiment. From this we hence conclude 
that the observed effects originate from the wave nature of the 
BEC and not from its many-body character. Nevertheless, ac- 
cording to Fig. S2, the mean-field potential stemming from 
the atomic interactions represents another important control- 
lable parameter that can be used to engineer the response of 
the system to the environmental action. 



